arXiv:1503.04185vl [cond-mat.mes-hall] 13 Mar 2015 


Minimal conductivity, topological Berry winding and duality in three-band semimetals 


Thibaud Louvet, Pierre Delplace, Andrei A. Fedorenko, and David Carpentier 

Laboratoire de Physique, Ecole Normale Superieure de Lyon, J^l allee dPtalie, 69007 Lyon, France 

(Dated: March 13, 2015) 


The physics of massless relativistic quantum 
particles has recently arisen in the electronic 
properties of solids following the discovery of 
graphene. Around the accidental crossing of two 
energy bands, the electronic excitations are de¬ 
scribed by a Weyl equation initially derived for 
ultra-relativistic particles. Similar three and four 
band semimetals have recently been discovered in 
two and three dimensions. Among the remark¬ 
able features of graphene are the characteriza¬ 
tion of the band crossings by a topological Berry 
winding, leading to an anomalous quantum Hall 
effect, and a finite minimal conductivity at the 
band crossing while the electronic density van¬ 
ishes. Here we show that these two properties 
are intimately related: this result paves the way 
to a direct measure of the topological nature of a 
semi-metal. By considering three band semimet¬ 
als with a flat band in two dimensions, we find 
that only few of them support a topological Berry 
phase. The same semimetals are the only ones 
displaying a non vanishing minimal conductivity 
at the band crossing. The existence of both a 
minimal conductivity and a topological robust¬ 
ness originates from properties of the underlying 
lattice, which are encoded not by a symmetry of 
their Bloch Hamiltonian, but by a duality. 

Electronic excitations which satisfy ultra relativistic 
quantum physics emerge around the band crossing of 
semimetals. While graphene constitutes the canonical 
example of such a phase, the discovery of topological in¬ 
sulators opened the route to other realizations [1]: such 
semimetals exist at the transitions between topological 
and trivial insulators. Theoretical proposals were initi¬ 
ated by the identification of crystalline symmetries stabi¬ 
lizing these critical phases both in two and three dimen¬ 
sions [2]. This idea turned out to be particularly suc¬ 
cessful: recently, a Weyl semimetal corresponding to two 
band crossings in three spatial dimensions was discov¬ 
ered in TaAs m- The existence of a three band cross¬ 
ing semimetal was proposed in a two dimensional carbon 
allotrope SG-lOb [5] and square M 0 S 2 sheet [6], and dis¬ 
covered in three dimensions in HgCdTe [7]. Stable four 
band crossing semimetals in three dimensions, denoted 
Dirac phases, were proposed theoretically and experi¬ 
mentally identified in NasBi [4l[8l|9] and in Cd 3 As 2 [T(H - 
[12], and predicted in /3-cristobalite Bi 02 [2] and distorted 
spinels m- 

In the case of graphene, the relativistic nature of the 
electronic excitations translates into remarkable trans¬ 
port properties including an anomalous half-integer Hall 


effect with half integer plateaus and a non vanishing min¬ 
imal conductivity exactly at the band crossing [HHizi. 
The anomalous half-integer Hall effect is related to the 
topological properties of the band crossing: when wind¬ 
ing around the crossing point, an electron picks up a 
quantized so-called Berry phase. This Berry phase is 
at the origin of the half-integer nature of the quantum 
Hall effect in graphene [m ng, but also characterizes 
the topological property of the semi-metallic phase en¬ 
coding its robustness towards gap opening perturbations 
[ 1 ]. On the other hand the minimal conductivity at the 
band crossing was associated to the Zitterbewegung of 
Dirac particles, an intrinsic agitation characteristic of ul¬ 
tra relativistic particles which leads to diffusive motion 
even in perfectly clean samples muni. 

In this article we consider these characteristic signa¬ 
tures of semimetals beyond the two-band crossing sit¬ 
uation of graphene. When a third locally flat band is 
present at the crossing, none of these properties is neces¬ 
sarily enforced. Indeed, we find that the exact same mod¬ 
els possess both a non vanishing minimal conductivity 
and a topological Berry winding at the crossing. Hence 
these properties are not hallmarks of a relativistic energy 
spectrum: they encode phase properties of the electronic 
wave functions. In the case of graphene, this phase wind¬ 
ing originates from the hexagonal lattice of carbon atoms. 
We will show that the lattice properties at the origin of 
these remarkable transport signatures are encoded not 
by a standard symmetry constraints of Bloch Hamilto¬ 
nians, but by a duality transformation. We identify two 
duality classes, corresponding to the two signatures: the 
existence or not of both a non-vanishing minimal con¬ 
ductivity and a topological Berry winding. 

Chirality or sublattice symmetry. We consider 
three-band semimetals in two dimensions characterized 
by two linearly crossing energy bands n = ± and a third 
locally flat band n = 0, represented in Fig. This 
corresponds to the situation of HgCdTe in three dimen¬ 
sions [7|, or the carbon allotrope SG-IOb [5] and square 
M 0 S 2 sheet [6| in two dimensions. The energy spectrum 
E{k) as a function of the momentum k exhibits the sym¬ 
metry E{k) —E{k) at least locally around the cross¬ 
ing point. This spectrum symmetry naturally originates 
from a chiral symmetry of the corresponding (low en¬ 
ergy) Hamiltonian. This symmetry is represented by a 
unitary operator C that anticommutes with the Hamil¬ 
tonian: PL = —CPLC. 

An explicit chiral operator can be defined when consid¬ 
ering the pedagogical examples of tight-binding Hamil¬ 
tonians defined on lattices. Similarly to graphene, only 


2 



FIG. 1 : Energy spectrum of a three-band chiral 
semimetal which shows a symmetry E{k) —E{k). 
Such a spectrum consists of two linear energy bands 
n = ± and a flat band n = 0. 


nearest neighbor couplings can be kept when focusing 
around the band crossing points. In this case chiral sym¬ 
metry corresponds to a sub-lattice symmetry: couplings 
are only present between the two sub-lattices A and B of 
a bipartite lattice. This is the case of the nearest neigh¬ 
bor description of graphene on the honeycomb lattice. In 
the case we consider in this article, three bands crossing 
implies the existence of three orbitals distributed on three 
Bravais lattices Ai , A2 and B of same geometry, as shown 
on Fig. Chiral symmetry is satisfied if the only cou¬ 
plings tiA2 relevant at low energy are between orbitals 
on the B and the ^1,^2 lattices whereas Ai and ^42 stay 
uncoupled. The corresponding Bloch Hamiltonian in the 
orbital basis (Ai, A2, B) is written 

^ / 0 0 \ 

H{ti,t2;k)= I 0 0 t2f2{k) ■ (1) 

Khftik) t2fUk) 0 / 


Such a Hamiltonian anti-commutes with a chirality oper¬ 
ator C = diag(l, 1 ,- 1 ). The complex functions fj{k) = 

\fj{k)\e'^^^^^^ encode the geometry of the lattice of cou¬ 
plings. Their amplitudes determine the spectrum of 
the semimetal: Eo(^) = 0 ^E±{k) = ±(t^|/i(^)p + 
^2|/2(^)P)^- a three band crossing occurs when fi and 
/2 vanish simultaneously at a point K in the Brillouin 
zone. In this article, we focus on properties which depend 
on the phases (l)j{k) and which are thus independent of 
the spectrum provided a band crossing occurs. 

Quite generally, we will consider chiral symmetric 
Bloch Hamiltonians describing a three-band semimetal 
with band crossing at point and whose linear expan¬ 


sion around the crossing takes the form 
H{K + = 

0 0 EiiPx + Ei2qy 

0 0 E 2 iqx + A22% 

^llQx + ^12Qy ^2 iQx + ^22% 0 

( 2 ) 

Such a Hamiltonian is entirely parametrized by a ma¬ 
trix A = {A^j} of complex coefficients. The phase of 
the coefficients A^j encodes the geometry of the under¬ 
lying lattice. The corresponding geometrical constraint 
on these phases must be independent of the amplitude 
of couplings between the orbitals. Hence it cannot re¬ 
sult in a symmetry of the Hamiltonian: we show that 
it corresponds to a duality in a manner analogous to the 
Kramers-Wannier duality of statistical mechanics models 
which relates models with different coupling amplitudes 
on different lattices |2U| . 

Lattice geometry and duality constraints. The 

geometry of the lattice can be described by two vectors 
ei and 62 relating a vector of the lattice B to neighbor 
sites of the Ai and A2 lattices, as shown on figureThe 
duality transformation V exchanges the lattices Ai and 
A2, or equivalently the vectors ei and 62, while simulta¬ 
neously exchanging the couplings between B and Ai or¬ 
bitals with couplings between B and A2 orbitals. Quite 
generally, this duality which is an involutive transforma¬ 
tion i.e. = I, relates a Hamiltonian H on a lattice jC 
to a Hamiltonian H on a different lattice C. However, 
on symmetric lattices where initial and dual lattices £, jC 
are related by a geometrical transformation IZ^ this du¬ 
ality translates into constraints on Hamiltonians defined 
on the same lattice (or same Hilbert space). In this case, 
and focusing for simplicity on nearest neighbor Hamilto¬ 
nians, the duality transformation can be recast into the 
form 



{DU)H{t 2 ,h-,TZk){DU)-^ = H{ti,t 2 -,k) ( 3 ) 

where 1 / is a unitary operator, IZ is the symmetry relating 
initial and dual lattices and D the operator swapping Ai 
and A2 orbitals: 

/O 1 0 \ 

D= 1 0 0 . ( 4 ) 

VO 0 Ij 

A very special case, which we call the duality class Vi 
corresponds to the situation where two orbitals lie on the 
same geometrical lattice, i.e. when ei = 62 7^ 0 . In this 
class, the duality transformation simplifies and U and IZ 
reduce to the identity. In this case, Bloch Hamiltonians 
encode the geometrical properties of a bipartite lattice, 
whereas in the other case, which we denote the duality 
class Pii, the underlying lattices are either Bravais lat¬ 
tices or possess three distinct sublattices. This duality 
restricts the form of the chiral tight-binding Hamiltonian 
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FIG. 2 : Triangular Bravais lattices with three orbitals (Ai, A2 and B) per unit cell. The position of the Ai orbitals 
is chosen arbitrary (a), in the center of the BAi hexagons (b), or at the same location than the A2 orbitals (c) The 
duality transformation V exchanges the location of the orbitals Ai and A2 and the hopping amplitudes symbolized 
by full/dashed lines. The original T3 lattice is recovered after an inversion in (b) whereas the dual and the original 

lattices coincide for the H3 lattice in (c). 


Q: in the class Vi we have fi(k) = f2{k) while a much 
weaker constrain fi{k) = f2{Bk) holds in class Vu for 
symmetric lattices. Generalizing this constrain to a local 
Bloch Hamiltonian (§ around a three band crossing, the 
duality in class Vi implies the condition 


A21 _ A22 
All Ai2 


( 5 ) 


while generically it only relates Hamiltonian at different 
crossing points in class Vu. In the following we show 
that Hamiltonians belonging to the duality class Vi de¬ 
scribe the only three-band semimetals possessing quan¬ 
tized topological Berry windings which are also those 
whose conductivity does not vanish at the band cross¬ 
ing and display a pseudo-diffusive regime. 

To illustrate this relation, let us discuss two nearest- 
neighbor lattice models of semi-metals belonging to 
both classes. A natural model in class Pj, inspired 
by graphene, corresponds to a honeycomb lattice with 
two orbitals on one of the two sub-lattices, shown on 
Fig. (c). The three bands of this model, which we 
call H3, cross at points K and K' of the Brillouin zone. 
Around those points, the Bloch Hamiltonian takes the 


form Q with a matrix of coefficients 



ti -iti \ 

t2 -it2 ) ’ 


(6) 


with a being the honeycomb lattice spacing, while the 
characteristic energy s cale of n earest neighbor couplings 
is encoded into t = This model satisfies the 

condition ^ with A = —i and belongs to the duality class 
Vi. The T3 model [ 21 ] is again defined on a honeycomb 
lattice but with additional orbitals A2 at the center of 
each hexagon as shown on Fig. (b). These orbitals 
are coupled by an amplitude t2 only to the B sub-lattice 
of the honeycomb lattice, while Ai and B orbitals on 
the honeycomb lattices are coupled by ti. This model 
belongs to the duality class Pn, with the inversion IZk = 
—k relating initial and dual lattices. Indeed, the Bloch 
Hamiltonian linearized around the band touching point 
K is written in the form ([^ with 


At 3 


3 a / h —iti \ 

Y\t2 it2 ) ' 


( 7 ) 


which does not fulfill the condition (§. Note that when 
ti = ^2, this linearized Hamiltonian can be written in 
the form H'K{q) = HvfS • where Sx, Sy and Sz = 
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diag(l, —1,0) satisfy the spin-1 algebra [ 5 '^, Sj] = icij^Sk- 
Hence the T3 model realizes a continuous deformation of 
spin-1 massless fermions, which all belong to the duality 
class Vii. 

We now characterize both the topological properties 
of the band crossing, as well as the electronic transport 
properties around the crossing, which turn out to be as¬ 
sociated to the duality class of the semi-metal. 


Topological property of a band crossing. The 

topological properties of band crossings in two dimen¬ 
sional crystals can be described by a topological invari¬ 
ant associated to each band around the crossing point. 
These invariants are the Berry phases 7^ (see eq.(lO) of 
Methods) acquired by Bloch eigenstates upon wind¬ 


ing around the crossing point K. By definition this Berry 
winding is independent of the path winding around the 
crossing point: it is an homotopic invariant of a given 
Hamiltonian. It describes a topological property of the 
semimetal when it is robust against perturbations of the 
Hamiltonian which do not lift the band crossing. Such a 
robustness occurs when this Berry winding is quantized. 
This is indeed the case in graphene where = ± 1 . For 
the three-bands semi-metals, these Berry windings can be 
readily obtained from the diagonalisation of the Hamilto¬ 
nian i)- We find that these windings are topological only 
when the condition (§ is fulfilled: the only semimetals 
characterized by a topological Berry winding are those of 
the duality class Pj, with values 


7+ = 7_ = sgnImA , 70 = —2sgnImA . (8) 


These quantized values of the Berry windings are stable 
with respect to any perturbation compatible with the du¬ 
ality constraint, i.e. which does not break the geometry 
of the underlying lattice. Eigenstates of the H3 model ^ 
are characterized by Berry phases 7± = — 1 , 70 = 2 as 
shown in Fig.[^ These windings are in particular robust 
to variations of hopping amplitudes ^1,^2- For any chiral 
symmetric semimetal which does not belong to the Vi 
class, the Berry phase can take any real value and de¬ 
pends continuously on deformations of the Hamiltonian. 
This result is illustrated for the T3 model 0 in Fig. 1 ^ 
the Berry phases 7^ are shown to vary continuously upon 
variation of the ratio ^2/^1 of nearest neighbor couplings. 

Minimal conductivity at the band crossing. 

Transport measurements constitute a powerful tool to 
probe the physical properties in the vicinity of the Fermi 
energy. We will show that close to the band crossing 
electronic transport is related to the phase content of 
the Hamiltonian (the phases 0 i(/c) of the amplitudes in 
Q) and not to the spectrum. Remarkably, in graphene, 
when the Fermi level coincides with the twofold band 
crossing point, the conductivity of a clean sample was 
predicted to remain finite despite a vanishing density of 
states. This result was first derived by considering the 
conductivity of a narrow strip of graphene between two 
contact electrodes as shown on Fig.|^ Let us consider an 
analogous setup for a three band chiral semimetal, i.e. 


a finite sample of length L and width W. Confinement 
of the sample between the leads gives rise to zero-energy 
evanescent states. At the band crossing, the conductiv¬ 
ity depends entirely on the nature of these evanescent 
states. Prior to an explicit calculation of the conduc¬ 
tivity (see Methods), it is instructive to consider the cur¬ 
rent operator jx{k) = defined from the tight- 

binding Bloch Hamiltonian Q. Introducing the ampli¬ 
tudes (7/^^^,'0^2 5 V^b) of fho electronic wavefunction in 
the three sub-lattices the longitudinal current can be ex¬ 
pressed as jx{k) = 2 Re['ipB{k){dk^fi{k)^\.{k))] and is 
found to be proportional to the amplitude on the 
B sub-lattice [ 23 ]. This hints that electronic transport 
at the band crossing will occur provided the zero-energy 
evanescent modes have a non-vanishing component on 
the B sub-lattice. 

As expected from the previous qualitative argument, 
the existence of a finite minimal conductivity at the 
threefold band crossing point is uniquely determined by 
the non vanishing weight of the wave function on the hub 
lattice as we have checked using a Landauer description 
of transport (see Methods). From this component is 
found to satisfy 

A.(»-)v.b = 0. (9) 

Hence the condition of existence of a non-vanishing min¬ 
imal conductivity at the band crossing is simply det A = 
0 . This constraint exactly coincides with the duality con¬ 
straint 0 defining the class Vp. the only three band 
semi-metals with a non-vanishing minimal conductiv¬ 
ity are exactly those belonging to this duality class Pj. 
Moreover for the H3 lattice model, this minimal conduc¬ 
tivity corresponds exactly to the value = e^/(7rh) 

predicted for graphene. This result remains valid for any 
model in the dual class Vi with an isotropic dispersion re¬ 
lation. For more general three band semi-metals with an 
anisotropic dispersion relation, the minimal conductivity 
depends on the angle between the leads and the princi¬ 
pal axes of the dispersion relation. However, the deter¬ 
minant of the corresponding conductivity tensor remains 
constant given by its isotropic value det a = {e‘^/(irh))^. 
In contrast, any three-band chiral symmetric semimetal 
that does not belong to the Vi class possesses a vanishing 
conductivity = 0 in every directions. Beyond the 

minimal conductivity, the fluctuations of this conductiv¬ 
ity can also be considered : their amplitude is encoded 
in the ratio between the shot noise power and the av¬ 
eraged current, the so-called Fano factor. This factor 
F takes a constant value F = 1/3 within the duality 
class Pi, a value already encountered in graphene m 
and characteristic of diffusive metals [ 24 ]. Such a result 
demonstrates that for all semi-metals in this class the 
transport through narrow perfectly clean junctions dis¬ 
plays the characteristic features of diffusive metals. 

We have evaluated the conductivity of different lattice 
models in the geometry of Fig. We compare the an¬ 
alytical results to a numerical Landauer approach (see 
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FIG. 3 : Momentum dependance of the phase associated with the Berry connection for the three energy levels 
n = —, 0, + for nearest neighbors models on the H3 and on the T3 model for three different values of the ratio of 
couplings between sub-lattices. A vortex is associated with a tt Berry phase. For the H3 lattice model, the Berry 
phase 7^ (10) along any closed loop winding around the origin is quantized (in units of tt) and signals a topological 
property of the band crossing points. This is not the case for the T3 lattice model, where this Berry phase 
continuously decreases and vanishes as a function of the couplings. 



FIG. 4 : The setup for two-terminal transport 
measurements: a 2 D sample of size IF x L with leads 
attached on opposite sides. 


Methods) to check for possible inter-crossing point ef¬ 
fects, neglected in the analytical approach. A perfect 
agreement is found between both approaches. The re¬ 
sults for the H3 model are shown as a function of the 


Fermi-energy, or gate potential Vg on Fig. 5 a the con¬ 
ductivity exhibits a plateau around the band crossing 
point = 0 , corresponding to cr = /{nh). The results 

of a similar study for the T3 model are also shown on 
Fig. 5 a and display a collapse of the conductivity around 
the band crossing point = 0 for three different val- 


ues of the couplings between orbitals (Inset). Fig. 5 b 


displays the analytical results for the dependance of the 
Fano factor on the gate potential Vg. We show that for 
the T3 model F(Vg = 0 ) = 1 , whereas for the H3 model 
the Fano factor reaches the value F = ^ at the band 
crossing, characteristic of a disordered metal as expected 
for class Pj. Finally, let us mention that for graphene the 
result of the Landauer formula for a narrow junction can 
be recovered for a long junction by an approach based on 
the Kubo formula m- While this equivalence remains 
valid for three band semimetals in the dual class Pj, it 
does not hold beyond it: we found that the Kubo con¬ 
ductivity for the T3 model diverges at the band crossing, 
as opposed to the vanishing Landauer conductivity for a 
narrow junction, in agreement with a previous result in 
the disordered limit 1^. 


As follows from our results, the occurrence of a trans¬ 
port regime with a non vanishing conductivity at the 
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FIG. 5 : Conductivity (a) and Fano factor (b) as a function of a gate potential Vg applied to the sample of size 
L = 100 , W = 300 in units of lattice spacing a. Fermi velocity is given by vp = 3 at/{ 2 h). Results for nearest 
neighbor lattices models on the H3 and T3 lattices are shown for various ratio ^2/^1 of couplings between the 

different lattices. 


crossing is not a generic property of linear dispersion 
relations near this crossing, nor a hallmark of relativis¬ 
tic physics of the associated electronic excitations such 
as the Zitterbewegung. It is indeed intimately related 
to the existence of a topological robustness of this band 
crossing, originating from lattice-properties encoded into 
the class Vi condition. This result opens the route to 
a direct probe of topological properties of semimetals 
through transport measurements around the band cross¬ 
ing. In graphene, a quantitative measurement of the con¬ 
ductivity exactly at the band crossing is hampered by 
fluctuations of the chemical potentials induced by a back 
gate. However the recent discovery of three dimensional 
semimetals changes the perspective: most of these new 
semi-metallic materials are stoichiometric and the Fermi- 
energy is expected to reside at the band crossing. In that 
situation the transport will be entirely dominated by the 
physics at the band crossing and will directly probe the 
topological properties associated with the band crossing 
points. 

Let us put in perspective the association between min¬ 
imal conductivity and topological properties raised by 
our results on three band semimetals in two dimen¬ 
sions. In two dimensions, we can extend our analy¬ 
sis for a band crossing described by the simple Hamil¬ 
tonian H{^ = SxQx + ^yQy^ with Sx^Sy satisfying a 
spin- 5 ' algebra. Again, we find that the Berry topo¬ 
logical winding m and the minimal conductivity are 
correlated. In particular for an integer spin 5 , corre¬ 
sponding to crossing of an odd number of bands, we find 
that the conductivity vanishes at the crossing, as does 
the Berry phase m, while both are finite for half in¬ 
teger spins. In three dimensions, a two band crossing 
is denoted a Weyl semimetal and is characterized by a 


topological Chern number instead of a Berry phase. In 
that case, the conductance at the nodal point remains 
finite and scales as G3d = ln2 ^/( 2 'kK){W /instead 
of G2d = I{ 7 rh){WlL) in two dimensions [ 28 ]. Sim¬ 
ilarly the value of the Fano factor is slightly modified 
from F2d = 1/3 to 0 . 57 . The existence of a fi¬ 

nite minimal conductance scaling as {W/L)^ appears to 
be a robust property associated with the topological two 
band crossing as its existence is neither modified by a 
weak disorder [ 55 ], nor by an anisotropic deformation of 
the cone or a tilt of this cone which breaks the local chiral 
symmetry m- 

METHODS 

Topological characterization. The topological 
characterization of a band crossing point is done by cal¬ 
culating the Berry phase associated to each eigenvec¬ 
tor |T^) upon winding anticlockwise around the crossing 
point: 

= . ( 10 ) 

Conductance and Fano factor. The conductance 
of a narrow sample is conveniently calculated from the 
set of the transmission probabilities of the conduc¬ 
tion channel labeled by n through the Landauer for¬ 
mula G = ^ Tn . The longitudinal conductivity a 
is related to this conductance as a = LW~^G. The 
Fano factor is related to these transmission coefficients 
as F = T^(l — T^)/ explicit calculation 

of the transmission probabilities requires solving the 
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Schrodinger equation piecewise and matching the solu¬ 
tions at the boundaries of the sample. 

Numerical calculations of transport were performed us¬ 
ing the Kwant code m, based on a Green function re¬ 
cursive technique to evaluate the transmission amplitude 
across a sample. Typical samples of dimensions L = 100 , 
W = 300 in lattice units were used, using a large po¬ 


tential in the electrodes eVooW/{vFh) = 45 , with Fermi 
velocity vf = 3 at/{ 2 h). 
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